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Abstract. A new Chemo-Dynamical Smoothed Particle 
Hydrodynamic (CD-SPH) code is presented. The disk 
galaxy is described as a multi-fragmented gas and star sys- 
tem, embedded in a cold dark matter halo with a rigid po- 
tential field. The star formation (SF) process, SNII, SNIa 
and PN events, and the chemical enrichment of gas, have 
all been considered within the framework of the standard 
SPH model, which we use to describe the dynamical and 
chemical evolution of triaxial disk-like galaxies. It is found 
that such approach provides a realistic description of the 
process of formation, chemical and dynamical evolution of 
disk galaxies over a cosmological timescale. 

Key words: chemical and dynamical evolution of disk 
galaxies - CD-SPH - star formation in SPH 



1. Introduction 

The dynamical and chemical evolution of galaxies is one of 
the most interesting and complex of problems. Naturally, 
galaxy formation is closely connected with the process of 
large-scale structure formation in the Universe. 

The main role in the scenario of large-scale struc- 
ture formation seems to be played by the dark matter. 
It is believed that the Universe was seeded at some early 
epoch with low density fluctuations of dark non-baryonic 
matter, and the evolving distribution of these dark ha- 
los provides the arena for galaxy formation. Galaxy for- 
mation itself involves collapse of baryons within poten- 
tial wells of dark halos flWhite fe Rees 1978|). The prop- 



erties of forming galaxies depend on the amount of bary- 
onic matter that can be accumulated in such halos and 
the efficiency of star formation. The observational sup- 
port for this galaxy formation scenario comes from the 
recent COBE detection of fluctuations in the microwave 
background (e.g. Bennett et al. 1993). 



The investigation of galaxy formation is a highly com- 
plex subject requiring many different approaches. The for- 
mation of self-gravitating inhomogeneities of protogalac- 
tic size, the ratio of baryonic and non-baryonic matter 
(Bardeen et al. 1986 



White & Silk 197£: Peebles 1993 



Dar 1995), the origin of the protogalaxy's initial angu- 



lar momentum (Voglis & Hiotelis 1989; Zurck et al. 1988 



Eisenstcin & Loeb 1995; Stcinmctz & Bartelmann 1995), 



and the protogalaxy's collapse and its subsequent evo- 
lution are all usually considered as separate problems. 
Recent advances in computer technology and numeri- 
cal methods have allowed detailed modeling of baryon 
matter dynamics in the universe dominated by colli- 
sionless dark matter and, therefore, the detailed gravi- 
tational and hydrodynamical description of galaxy for- 
mation and evolution. The most sophisticated models 
include radiative processes, star formation and super- 
nova feedback (e.g. [Katz 199$ [Stcinmctz fc Mullcr 1994; 



Fricdli fc Bcnz 1995) 
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The results of numerical simulations are essentially af- 
fected by the star formation algorithm incorporated into 
modeling techniques. The star formation and related pro- 
cesses are still not well understood on either small or large 
spatial scales, such that the star formation algorithm by 
which the gas material is converted into stars can only be 
based on simple theoretical assumptions or on empirical 
observations of nearby galaxies. The other most impor- 
tant effect of star formation on the global evolution of a 
galaxy is caused by a large amount of energy released in 
supernova explosions and stellar winds. 

Among numerous methods developed for the modeling 
of complex three dimensional hydrodynamic phenomena, 
Smoothed Par ticle Hydrodyna mics (SPH) is one of the 
most popular ( Monaghan 1992 ). Its Lagrangian nature al- 
lows easy combination with fast N-body algorithms, mak- 
ing it suitable for simultaneous desc ription of the comple x 
dynamics of a gas-stellar system ( Friedli fc Bcnz 1995 ). 
As an example of such a com bi nation, the TREE-SPH 
code ( |Hernquist fc Katz 1989| ; |Navarro fc White 1993b 
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was successfully applied to the detailed modeling of disk 
galaxy mergers ( Mihos fc Hcrnquist 199(f ) and of galaxy 
formation and evolution (Katz 1992). The second good ex- 
ample is a GRAPE-SPH code (|Stcinmctz fc Muller 1994; 



neighborhood of each particle within the 2 • hi remains 



Steinmctz & Muller 1995) which was successfully used to 



model the evolution of disk galaxy structure and kinemat- 
ics. 

In recent years, there have been excellent papers con- 
cerning the complex SPH modeling of galaxy formation 
and evolution (Raiteri et al. 1996; Carraro et al. 1997). 
Our code, proposes new "energetic" criteria for SF, and 
suggest a more realistic account of returned chemically 
enriched gas fraction via SNII, SNIa and PN events. 

The simplicity and numerical efficiency of the SPH 
method were the main reasons why we chose this tech- 
nique for the modeling of the evolution of complex, multi- 
fragmented triaxial protogalactic systems. We used our 
own modification of the hybrid N-body/SPH method 
(Berczik & Kravchuk 1996; |Berczik 199§| ), which we call 
the chemo-dynamical SPH (CD-SPH) code. 

The "dark matter" and "stars" are included in the 
standard SPH algorithm as the N-body collisionlcss sys- 
tem of particles, which can interact with the gas com- 
ponent only through gravitation ( Katz 1992] ). The star 
formation process and supernova explosions are also in- 
cluded in the scheme as proposed by Raiteri et al. (1996), 
but with our own modifications. 



2. The CD-SPH code 

2.1. The SPH code 

Continuous hydrodynamic fields in SPH are described by 
the interpolation functions constructed from the known 
values of these functions at randomly positioned parti- 
cles ( Monaghan 1992 ). Following Monaghan & Lattanzio 
(1985) we use for the kernel function Wij the spline ex- 
pression in the form of: 



f 1 - §tt?- + fufj, if < Uij < 1, 



1 



13 nh 3 | 4 


V 

Here u 



if 1 < Uij < 2, 



otherwise . 



(1) 



' 1 1 — 1 ij / h. 

To achieve the same level of accuracy for all points in 
the fluid, it is necessary to use a spatially variable smooth- 
ing length. In this case each particle has its individual 
value of h. Following Hernquist & Katz (1989), we write: 



N 1 



(2) 



In our calculations the values of hi were determined 
from the condition that the number of particles Nb in the 



constant (Mihos & Hcrnquist 1996). The value of Nb is 
chosen such that a certain fraction of the total number 
of "gas" particles N affects the local flow characteristics 



(Hiotelis & Voglis 1991). If the defined hi becomes smaller 
than the minimal smoothing length h m i n , we set the value 
hi = h min . For "dark matter" and "star" particles (with 
Plummer density profiles) we use, accordingly, the fixed 
gravitational smoothing lengths hd m and h star - 

If the density is computed according to Equation (||), 
then the continuity equation is satisfied automatically. 
Equations of motion for particle i are 



dt 

dvj 
~dt 



Pi 



(3) 



(4) 



where Pi is the pressure, $j is the self gravitational 
potential, <&1 xt is a gravitational potential of possible 
external halo and is an artificial viscosity term 

(Hiotelis et al. 1991). The energy equation has the form: 

dui Pi Ti~ A t 

-jr = V,v, + . (5) 

at pi p t 

Here Ui is the specific internal energy of particle i. The 
term (P, — hi)/ pi accounts for non adiabatic processes not 
associated with the artificial viscosity (in our calculations 
F, = 0). We present the radiative cooling in the form: 



A,; = hi(ui,pi) = h*(Ti) ■ n\ 



(6) 



where rii is the hydrogen number density and Ti the 
temperature. To follow its subsequent thermal behaviour 
in numerical simulations, we use an analytical approxi- 
mation of the standard cooling function A*(T) for an op- 
t ically thin primordial plasma in ionization e quilibrium 
( |Dalgarno fc McCray 1972| ; |Katz fc Gunn 199l|) . Its abso- 
lute cutoff temperature is set equal to 10 4 K. 

The equation of state must be added to close the sys- 
tem. 



Pi = Pi ■ (7 - 1) • Ui, 

where 7 = 5/3 is the adiabatic index. 



(7) 



2.2. Time integration 

To solve the system of E quations (p|), (p|) and (p[ ) we 
use the leapfrog integrator ( Hcrnquist fc Katz 1989[). The 



time step Sti for each particle depends on the particle's ac- 
celeration an and velocity v^, as well as on viscous forces. 
To define Sti we use the relation from Hiotelis & Voglis 
(1991), and adopt C ourant's number C n =0. 1. 

We carried out ( Berczik fc Kolesnik 1993| ) a large se- 
ries of test calculations to check that the code is correct, 
the conservation laws are obeyed and the hydrodynamic 
fields are represented adequately, all with good results. 
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2. 3. The star formation algorithm 



where: 



N e 



It is well known that star formation (SF) regions are as , , 

sociated with giant molecular complexes, especially with v c = ^ m j ' v j/ ^ 
regions that are approaching dynamical instability. The 
early phase of star formation docs not seem to crucially 
affect the dynamics of a galaxy. From the beginning of 
the collapse, such a system decouples from its surround- 
ings and evolves on a completely different timescale. When 
the chemically enriched gas content of the galaxy de- 
creases, the heating by winds and supernova explosions 
(Lcithcrer et al. 1992) begins to play an important role 
in the dynamics of the galaxy. The overall picture of 
star formation seems to be understood, but the detailed 
physics of star formation and accompanying processes, on 
cither sma ll or large scales, remains sketchy ( Larson 19(>9[ 



(11) 



For practical reasons, it is useful to define a critical 
temperature for SF onset in particle i as: 



(12) 



Then, if the temperature of the "gas" particle i, drops 
below the critical one, SF can proceed. 



T, < Tf 



(13) 



Silk 1987|) 



All the above stated as well as computer constrains 
cause the using of simplified numerical algorithms of de- 
scription of conversion of the gaseous material into stars, 
which are based on simple theoretical assumptions and/or 
on results of observations of nearby galaxies. 

To describe of the process of converting of gaseous ma- 
terial into stars we modify the standard SPH star f orma- 
tion algorithm flKatz 1992^ |Navarro fc White 1993|), tak- 



ing into account the presence of random motions in the 
gaseous environment and the time lag between the ini- 
tial development of suitable conditions for star forma- 
tion and star formation itself ( Berczik fe Kravchuk 199(| 
Berczik 1998). The first reasonable requirement incorpo- 
rated into this algorithm allows selecting "gas" particles 
that are potentially eligible to form stars. It states that in 
the separate " gas" particle the SF can start if the absolute 
value of the " gas" particle's gravitational energy exceeds 
the sum of its thermal energy and the energy of random 
motions: 



Ef \>E* h + Ef h . 



(8) 



Gravitational and thermal energies and the energy of 
random motions for the "gas" particle i in model simula- 
tion are defined as: 



Ef = -\.G.m\lK 
Ef = | • m, ■ cl 



(9) 



Ef = \ ■ m t ■ Avl 



where Ci = ^/3? • Tij /i is the isothermal sound speed 
of particle i. We set fi = 1.3 and define the random or 
"turbulent" square velocities near particle i as: 



Ne 



Avf =J2 m j- (vj - v c ) 2 / ' 



(10) 



We think that requirement (|8j), or in another form (|13[), 
is the only one needed. It seems reasonable that the chosen 
"gas" particle will produce stars only if the above condi- 
tion hold over the interval that exceeds its free - fall time 
tfj — \/3 ■ 7r/ (32 • G ■ p). This condition is based on the 
well known fact that, due to gravitational instability, all 
substructures of a collapsing system are formed on such a 
timescale. Using it, we exclude transient structures, that 
are destroyed by the tidal action of surrounding matter 
from consideration. 

We also define which "gas" particles remain cool, 
i.e. t coo i < tff. We rewrite this condition as presented 
in Navarro & White (1993): pi > p cr it- Here we use the 
value of p cr u — 0.03 cm~ 3 . 

When the collapsing particle i is defined, we create the 
new "star" particle with mass m star and update the "gas" 
particle rrii using these simple equations: 



m star = e • m,i, 



mi = (1 - e) • mi 



(14) 



Here e, defined as the global efficiency of star forma- 
tion, is the fraction of gas converted into stars according 
to the appropriate initial mass function (IMF). The typ- 
ical values for SF efficiency in our Galaxy on the scale 
of giant molecular clouds are in the range e « 0.01 -j- 0.4 



(Ducrr et al. 1982; |Wilking fc Lada 1983; ). But it is still a 
little known quantity. In numerical simulation the model 
parameter has to be checked by comparison of numerical 
simulation results with available observational data. Here 
we define e as: 



e = l-(E^+Ef)/\E 



(15) 



with the requirement that all excess mass of the gas 
component in a star-forming particle, which provides the 
inequality | Ef |> E\ h + Ef h , is transformed into the 
star component. In the code we set the absolute maxi- 
mum value of the mass of such a "star" particle = 
2.5 10 6 M© i.e. « 5% of the initial particle mass m^. 
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At the moment of stellar birth, the position and veloc- 
ities of new "star" particles are equal to those of parent 
"gas" particles. Thereafter these "star" particles interact 
with other " gas" , " star" or " dark matter" particles only by 
gravitation. The gravitational smoothing length for these 
(Plummer like) particles is set equal to h s tar- 

2-4- The thermal S Nil feed-back 

We try to include the events of SNII, SNIa and PN in the 
complex gasdynamic picture of galaxy evolution. But, for 
the thermal budget of the ISM, only SNII plays the main 
role. Following Katz (1992), we assume that the explo- 
sion energy is converted totally into thermal energy. The 
stellar wind action seems not to be essential in the en- 



ergy budget ( [Ferriere 1995[ ) . The total energy released by 
SNII explosions (10 44 J per SNII) within a "star" particle 
is calculated at each time step and distributed uniformly 
between the surrounding (i.e. Tij < h sta r) "gas" particles 
(Raiteri et al. 1996). 



'oo(Z) = 10.130 + 0.0755 -logZ- 0.0081 • 
ox(Z) = 4.4240 + 0.7939 ■ logZ + 0.1187 ■ 
a 2 (Z) = 1.2620 + 0.3385 -log Z + 0.0542- 



(logZ) 2 , 
(logZ) 2 , 
(logZ) 2 . 



(18) 



These relations are based on the calculations of 
the Padova group (Alongi et al. 1993; Bressan et al. 1993; 
Bcrtelli et al. 1994) and give a reasonable approximation 
to stellar lifetimes in the mass range from 0.6 Mq to 
120 M Q and metallicities Z = 7 • 10~ 5 0.03 (defined 
as a mass of all elements heavier than He). In our cal- 
culation following Raiteri et al. (1996), we assume that Z 
scales with the oxygen abundance as Z/Z Q = 16 O/ 16 O . 
For those metallicities exceeding available data we take 
the value corresponding to the extremes. 

We can define the number of SNII explosions inside a 
given "star" particle during the time from t to t + At using 
a simple equation: 



2. 5. The chemical enrichment of gas 

Every "star" particle in our SF scheme represents a sep- 
arate, gravitationally closed star formation macro region 
(like a globular cluster). The "star" particle is character- 
ized by its own time of birth tb eg SF which is set equal 
to the moment of particle formation. After the formation, 
these particles return the chemically enriched gas into sur- 
rounding "gas" particles due to SNII, SNIa and PN events. 
For the description of this process we use the approxima- 
tion proposed by Raiteri et al. (1996). We consider only 
the production of 16 and 56 Fe, and try to describe the 
full galactic time evolution of these elements, from the 
beginning to present time (i.e. t evo i w 13.0 Gyr). 

With the multi-power IMF law suggested by Kroupa 
et al. (1993), the distribution of stellar masses within a 
"star" particle of mass m star is then: 



*(m) = m star ■ A • < 



777 



-2.7 



m- 1 - 3 , if 0.1 < m < 0.5, 

! > if 0.5 < 777 < 1.0, (16) 

if 1.0 < 777, 



where m is the star mass in solar units. With adopted 
lower (mi ow = 0.1 Mq) and upper (m upp — 100 Mq) 
limits of the IMF, the normalization constant A w 0.31. 

For the definition of stellar lifetimes we use the equa- 
tion (Raiteri et al. 1996): 

\ogt dead = oo(Z) - oi(Z) • log 777 + a 2 (Z) • (logm) 2 , (17) 

where tdead is expressed in years, m is in solar units, 
and coefficients arc defined as: 



AA SN „ = 



^(m)dm, 



(19) 



m dcad (t+At) 



where mdead(t) and mdead(t + At) are masses of stars 
that end their lifetimes at the beginning and at the end of 
the respective time step. We assume that all stars with 
masses between 8 Mq and 100 Mq produce SNII, for 
which we use the yields from Woosley & Weaver (1995). 
The approximation formulae from Raiteri et al. (1996) de- 
fines the total ejected mass by one SNII - m*°*, as well as 



the ejected mass of iron ■ 



Fc 



and oxygen 



tion of stellar mass (in solar units). 



to?, as a func- 

C J 



<m\f 



7.682- 10" 1 -to 



1.056 



m\° = 2.802- 10~ 4 -m 1864 , 



TO o. = 4.586 -10- 4 -777 2 - 721 . 



(20) 



To take into account PN events inside the "star" par- 
ticle we use the equation, as for (|l^): 



AAp N = J $(m)dm. 

m dBad ,(t+At) 



(21) 



Following van den Hoek & Groenewegen (1997), Sam- 
land (1997) and Samland et al. (1997), we assume that 
all stars with masses between 1 Mq and 8 Mq pro- 
duce PN. We define the average ejected masses (in 



solar units) of one PN event as 


Rcnzini & Voli 1981 


van den Hock & Groenewegen 1997 


: 
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mlf = 1.63, 



m\] = 0.00, 



(22) 



0.(10. 



The method described in Raiteri et al. (1996) and pro- 
posed in Greggio & Renzini (1983) and Matteuchi & Greg- 
gio (1986) is used to account for SNIa. In simulations, the 
number of SNIa exploding inside a selected " star" particle 
during each time step is given by: 



ATVsNIa = J $2(^2)^2- 

rn dead (t+At) 



(23) 




■ ' 



Fig. 1. The number of SNII, SNIa and PN events 



The quantity ^2(^2) represents the initial mass func- 
tion of the secondary component and includes the distri- 
bution function of the secondary's mass relative to the 
total mass of the binary system tub, 



'"'sup 

* 2 (m 2 ) =m star -A 2 - [ (— J^m/V 
J m B 



(24) 



where m m j = max(2 • to 2 , 3 Mq) and m sup = rn 2 + 
8 M Q . Following van den Berg & McClure (1994) the value 
of normalization constant we set, equal to A2 = 0.16 • A. 

The total 
ejected mass (in solar units) is (Thielemann et al. 1986; 
Nomoto et al. 1984): 



..tot 

1% = 0.13. 



0.63, 



(25) 



In summary, a new "star" particle (with metallicity 
Z = 10~ 4 ) with mass 10 4 Mq during the total time of 
evolution t evo i produces: 



r Amg NII 


= 477, 


Am£ N 


= 2164, 


Am SNIa 


= 4.14, 


Am SNII 


= 159, 


Am^ 


= 721.3, 


AmfSia 


= 1.38, 


Am SNII 


= 3.5, 


Am{^ 


= 0.000, 


A ™SNIa 


= 5.35, 


, Am SNII 


= 119, 


Amg N 


= 0.000, 


ATO SNIa 


= 1.10. 



(26) 




AiVgNn « 52.5, AiV PN « 1770.0, AiV SNIa w 8.4 



Fig. 2. The returned mass of 56 Fe 



Fig. 0. presents the number of SNII, SNIa and PN 
events for this "star" particle. 

In Fig. |2. and Fig. ||. we present the returned masses 
of 56 Fe and 16 0. We can estimate the total masses (H, He, 
56 Fe, 16 0) (in solar masses) returned to the surrounding 
"gas" particles due to these processes as: 



2.6. The cold dark matter halo 

In the literature we have found some profiles, some- 
times controversial, f or the galact i c Cold Dark Mat- 
ter Haloes (CDMH) ( |Burkert 1995| ; |Navarro 19981 ). For 



resolved structures of CDMH: Phalo{r) 
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Fig. 3. The returned mass of ie O 

(Moore et al. 1997). The structure of CDMH high 
resolution N-body simulations can be described by: 
Phaio(r) ~ r^ 1 (Navarro et al. 1996; Navarro et al. 1997). 
Finally, in Kravtsov et al. (1997), we find that the cores 
of DM dominated galaxies may have a central profile: 
Phaio{r) ~ r-° 2 . 

In our calculations and as a first approximation, it 
is assumed that the model galaxy halo contains the 
CDMH component with Plummer - type density profiles 
( [Doupholc fc Colin 199E| ): 



Phalo{r) 



halo 



b 5 

u halo 



(27) 



Therefore for the external force acting on the "gas" 
and "star" particles we can write: 



-G ■ 



halo 



(28) 



3. Results and discussion 

3.1. Initial conditions 

After testing our code demonstrated that simple assump- 
tions can lead to a reasonable model of a galaxy. The 
SPH calculations were carried out for N gas = 2109 "gas" 
particles. According to Navarro & White (1993) and Rai- 
teri et al. (1996), such a number seems adequate for a 
qualitatively correct description of the systems behaviour. 
Even this small a number of "gas" particles produces 
Nstar = 31631 "star" particles at the end of the calcu- 
lation. 

The value of the smoothing length hi was chosen to 
require that each "gas" particle had Nb — 21 neighbors 
within 2 ■ hi. Minimal /i m -m was set equal to 1 kpc, and 
the fixed gravitational smoothing length h star — 1 kpc 
was used for the "star" particles. Our results show that 



a value of Nb ~ 1% N gas provides qualitatively correct 
treatment of the system's large scale evolution. 

As the initial model (relevant for CDM - scenario) 
we took a constant-density homogeneous gaseous triaxial 
configuration (M gas = 10 11 M Q ) within the dark matter 
halo (M halo = 10 12 M Q ). We set A = 100 kpc, B = 75 
kpc and C — 50 kpc for semiaxes of system. Such triaxial 
configurations are reported in co smological simulations of 



the dark matter halo formation ( Eiscnstcin fc Loeb 1995 ; 
Frenk et al. 1988; Warren et al. 1992). Initially, the cen- 
ters of all particles were placed on a homogeneous grid 
inside this triaxial configuration. We set the smoothing 
parameter of CDMH: 6/ la / = 25 kpc. These values of 
Mhaio and btaio are typical for CDMH in disk galaxi es 
(Navarro et al. 1996; Navarro et al. 1997; [Burkert 1995| ). 

The gas component was assumed to be cold initially, 
To = 10 4 K. As we see in our calculations, the influence 
of random motions essentially reduces the dependence of 
model parameters on the adopted temperature cutoff and, 
therefore, on the adopted form of the cooling function it- 
self. 

The gas was assumed to be involved in the Hubble 
flow (Hq = 65 km/s/Mpc) and the solid - body rotation 
around z - axis. We added small random velocity compo- 
nents (A | v |= 10 km/s) to account for the random mo- 
tions of fragments. The initial velocity field was defined 



v(x, y, z) = [fl(x, y 1 z)xr]+H -r + Av(x, y, z), (29) 

where £l(x,y, z) is the angular velocity of an initially 
rigidly rotating system. 

The spin parameter in our numerical simulations is 
A « 0.08, defined in Peebles (1969) as: 



G ■ {M gas + M halo )^ 



(30) 



Lo is the total initial angular momentum and Eq T is 
the total initial gravitational energy of a protogalaxy. It 
is to be noted that for a system in which angular mo- 
mentum is acquired through the tidal torque of the sur- 
rounding matter, the standard spin parameter d oes not 
exceed A « 0.11 ( Steinmetz k, Bartelmann 1995 ). More- 
over, its typical values range between A w 0.07;^;^ 4 , 
e.g. 0.02 < A < 0.11. 

3.2. Dynamical model 

In Fig. |. we present the "XY", "XZ" and "YZ" distribu- 
tions of "gas" particles at the final time step (t evo i « 13.0 
Gyr). The box size is 50 kpc. In Fig. 0. we present the 
distributions of "star" particles. The "star" distributions 
have dimensions typical of a disk galaxy. The radial ex- 
tension is approximately 25 — 30 kpc. The disk height is 
around 1 — 2 kpc. In the center the "bar -like" structure 
is developed as a result of strong initial triaxial structure 
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whole in the plane of the disk we can see the " spiral - like" 
distribution of particles, with extended arm filaments. The 
"gas" particles are located within central 5 — 10 kpc. 

Except for the central region (< 2 kpc), the gas dis- 
tribution has a exponential form with radial scale length 
w 2.8 kpc. The column density distributions of gas a gas {r) 
and stars er* (r) are presented in Fig. |[ The total column 
density is defined as: <J to t{ r ) = c gas (r) + o*(r). The to- 
tal column density distribution atot{r) is well approxi- 
mated (in the interval from 5 kpc to 15 kpc) with an 
exponential profile characterized by a ss 3.5 kpc radial 
scale length. This value is very close to one a reported 
recently (3.5 kpc) for the radial scale length of the to- 
tal disk mass surface density distribution obtained for our 
Galaxy (Mera et al. 1998b). The value of a tot w 55 Mq 
pc -2 near the location of the Sun (r « 9 kpc) is close to a 
recent determination of the total density 52 ±13 Mq pc~ 2 
(Mera et al. 1998a). 

Fig. ^. shows both the rotational velocity distribu- 
tion of gas V ro t{r) resulting from the modeled disk 
galaxy calculation and the rotational curve for our Galaxy 



(Vallce 1994), both of which are very close. 

The gaseous radial V ra d{r) and normal V z (r) velocity 
distributions are in Fig. M. and Fig. ^|. The radial veloc- 
ity dispersion has a maximum value « 60 km/s in the 
center, a high value mainly caused by the central strong 
bar structure. Near the Sun this dispersion drops down 
to Ri 20 km/s. Such radial dispersion is reported in the 
kinematic study of the stellar motions in the solar neigh- 
borhood ( |Bicnaymc 199§| ), while the normal dispersion is 
near ps 20 km/s in the whole disk. This value also coin- 
cides with the vertical dispersion velocity near the Sun 



( [Bicnaymc 1998[ ). 

We present the temperature distribution of gas T(r) in 
Fig. [n]. As seen the distribution of T(r) has a very large 
scatter from 10 4 K to 10 6 K. In our calculation we set the 
cutoff temperature for the cooling function at 10 4 K, the 
gas can't cool to lower temperatures. 

The modeled process of SNII explosions injects to a 
great amount of thermal energy into the gas and gener- 
ates a very large temperature scatter, also typical of our 
Galaxy's ISM. At each point even with crude numerical 
approximations a good fit can be reached for all dynami- 
cal and thermal distributions of gas and stars in a typical 
disk galaxy like our Galaxy. 

3.3. Chemical characteristics 



Fig. O. shows the time evolution of the SFR in galaxy 
SFR(t) = dM*{t)/dt. Approximately 90% of gas is con- 
verted into stars at the end of calculation. The most in- 
tensive SF burst happened in the first « 1 Gyr, with a 
maximal SFR w 35 M Q /yr. After w 1.5 - 2 Gyr the SFR 
is decreases like an "exponential function" until it has a 
value « 1 Mq I yr at the end of the simulation. To check 
the SF and chemical enrichment algorithm in our SPH 




Fig. 6. (T g as(r), a*(r) and a to t(r) = o gaa {r) + cr»(r). The 
column density distribution in the disk of gas and stars in the 
final step 

300.0 



250.0 -_ **fj. : 



3 i oO, Q , r-°ot)b 



oo 



o - Vallee (1994) 



:1 



0.0 5.0 10.0 15.0 20.0 25.0 

r (kpc) 



Fig. 7. V ro t(r). The rotational velocity distribution of gas in 
the final step 



code, we use the chemical characteristics of the disk in 
the "solar" cylinder (8 kpc < r < 10 kpc). 

The age-metallicity relation of the "star" particles in 
the "solar" cylinder, [Fe/H](i), is shown in Fig. |l2L with 
observational data taken from Meusinger et al. (1991) and 
Edvardsson et al. (1993), while in Fig. [l3|. we presented 
the metallicity distribution of the "star" particles in the 
"solar" cylinder A*([Fe/fi]). The model data are scaled to 
the observed number of stars (Edvardsson et al. 1993). In 
Fig. O. each model point represents the separate "star" 
particle. The mass of each "star" particle is different (from 
~ 10 4 Mq up to ~ 2.5 10 6 Mq), because the star forma- 
tion efficiency - e is different in each star forming region. 
The model point is systematically higher than the obser- 
vations (especially near the t ss 5 Gyr), but if we also 
analyze the mass of each "star" particle we see that the 
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Fig. 4. The distribution of "gas" particles in the final step 



more massive particles systematically show lower metal- 
licity than the observations. If one divides the metallicity 
to equal zone and calculate the sum of the mass in each 
metallicity zone in Fig. we get the results in Fig. |l3|. 
and in this figure we see what the model mass distribution 
has shifted to the lower metallicities. 

The [O/Fe] vs. [Fe/H] distribution of the "star" par- 
ticles in the "solar" cylinder one found in Fig. [l4]. In this 
figure we also present the observational data from Ed- 
vardsson et al. (1993) and Tomkin et al. (1992). All these 
model distributions are in good agreement, not only with 
presented observational data, but also with other data col- 
lected from Portinari et al. (1997). 

The [O/H] radial distribution [0/H](r) is shown in 
Fig. [l5|. The approximation presented in the figure is ob- 
tained by a least - squares linear fit. At distances 5 kpc 
< r < 11 kpc the models radial abundance gradi- 
ent is —0.06 dex/kpc. In the literature we found differ- 
ent values of this gradient defined in obj ects of different 
types. From observations of HII regions flPcimbcrt 197§ 
Shaver et al. 1983) we obtained oxygen radial gradient 



0.07 dex/kpc. From observa tions of PN of different 



types ( Maciel fc Koppcn 1994 ) we obtained the values: 
-0.03 dex/kpc for PNI, -0.069 dex/kpc for PNII, -0.058 
dex/kpc for PNIII, -0.062 dex/kpc for PNIIa, and -0.057 
dex/kpc for PNIIb. All this agrees well with the oxygen 
radial gradient in our Galaxy. 

3.4. Conclusion 

This simple model provides a reasonable, self-consistent 
picture of the processes of galaxy formation and star for- 
mation in the galaxy. The dynamical and chemical evo- 
lution of the modeled disk-like galaxy is coincident with 
observations for our own Galaxy. The results of our model- 
ing give a good base for a wide use of the proposed SF and 
chemical enrichment algorithm in other SPH simulations. 
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